Depolarization of backscattered linearly polarized light 
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We formulate a quantitative description of backscattered linearly polarized 
light using an extended photon diffusion formalism taking explicitly into ac- 
count the scattering anisotropy parameter g of the medium. From diffusing 
wave spectroscopy measurements the characteristic depolarization length for 
linearly polarized light l v is deduced. We investigate the dependance of this 
length on the scattering anisotropy parameter g spanning an extended range 
from -f (backscattering) to 1 (forward scattering). Good agreement is found 
with Monte Carlo simulations of multiply scattered light. 
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1. Introduction 

Polarized light scattered many times in a random 
medium leaves the sample partially depolarized. Unfor- 
tunately, despite its importance in areas like biomedi- 
cal optical imaging, coherent backscattering or dynamic 
spectroscopjiSi 3 ^, the depolarization of light in a ran- 
dom medium is still not completely understood due to 
the complexity of vector wave multiple scattering (as 
compared to the much simpler problem of scalar wave 
propagation). Previous attempts have mainly focussed 
on isotropic (Rayleigh) scattering 5 , or on the depolariza- 
tion of circularly polarized lights^. Only few studies have 
discussed specifically the mechanism of depolarization of 
linearly polarized light in the case where the anisotropy 
parameter g =< cos (9 > is different from and further- 
more a detailed comparison with experiment has been 
lacking§iL§i2ii2iii. An accurate description for arbitrary 
scattering anisotropy is however crucial to analyze the 
information contained in backscattered light if progress 
is to be made in applications like remote sensing, photon 
correlation spectroscopy or optical imaging of biological 
tissuea 1 ! 13 ! 14 . 

In this paper we formulate a quantitative description of 
backscattered linearly polarized light using an extended 
photon diffusion formalism taking explicitly into account 
the scattering anisotropy parameter g. The details of 
our model are adjusted by comparison with Monte Carlo 
simulations of multiply scattered light. We show how the 
characteristic length of depolarization of incident linearly 
polarized light i v can be deduced from measurements of 
intensity fluctuations of light scattered from liquid turbid 
media via Diffusing Wave Spectroscopy (DWS). We can 
distinguish the following limiting situations for the trans- 



port of light and its polarization: isotropic scattering 
g ~ 0, forward-peaked scattering g ~ 1 and backward- 
peaked scattering g ~ —1. The situation of forward- 
peaked scattering g ~ 1 is typical of Mie scatteringi 5 . 
with large particles and of biological tissues, whereas the 
situation of g < has only been made possible experi- 
mentally recently by tuning the interaction of the light 
using mesostructured colloidal liquids^. Here we show 
that with our additional correction the simple photon 
diffusion picture successfully describes the distribution 
of path lengths and the DWS autocorrelation function in 
the backscattering geometry. 

2. Path lengths distribution for backscattering 

On length scales much larger than the transport mean 
free path I*, the transport of light in a turbid medium 
can be described by the diffusion approximation. This 
approximation is connected to the idea of treating the 
transport of photons as a random walk, characterized by 
a distribution of path lengths&i£>i& 1 £. An exact solution 
of the diffusion equation applied to light transport can 
be obtained using the method of images. This method 
takes into account the boundary conditions through two 
lengths which are both of the order of a transport mean 
free path: the extrapolation length z e where the flux of 
the flux of photons vanishes outside the sample and z p 
which is the location of the photon source (for a more 
detailed interpretation of z e , z p based on a random walk 
model see reference 4) . The method leads to 



P(s) = 



V3 



+ {z p + 2z e ), 



3 (zp+2; 
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(1) 

which obeys the normalization condition J Q P(s)ds = 1. 
Note that the path length is simply related to the number 
of scattering events n by s/l = n — 1, so that the path 
length is for single scattering. Here I is the scattering 
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anisotropic scattering corresponding to g ~ 0.806. 



f(9) = 3g/2 




Fig. 1. Normalized path length distribution P(s) x I for 
backscattered light from a semi-infinite medium. Sym- 
bols: Monte Carlo simulations. Lines: calculations based 
on Eq.EJwith f(g) = 3g/2 (solid) and f{g) = (dashed). 
Inset: f(g) obtained from Eq.^adjusted to fit the simu- 
lation results. Wavelength Ao = 532 nm, refractive index 
of the particle n p = 1.59 and of the solvent n s = 1.332. 
Non-reflecting boundary conditions were used. 



• 500 nm,l t = 4.8 1, g=0.S56 
O 400 nm, I = 3.6 I, g=0.806 
■ 80 nm, I =2.6 I, g=0.068 




Fig. 2. Depolarization of multiply scattered light. 
2d«(s) — 1 from theory (Eq.J3J solid lines) and simula- 
tion (symbols) . Excellent agreement if found for Rayleigh 
scatterers while for larger particles the agreement be- 
comes somewhat less good. 

The disagreement in the latter case is not surprising as 
the diffusion approximation is known to overestimate the 
contribution from the short paths of the distribution, the 
error becoming more and more severe as the anisotropy 
of scattering increases. One way to improve the distri- 
bution of paths length of Eq. ^ is by introducing a cut- 
off in the distribution as suggested by Mackintosh and 
Johni Here we extend their approach by taking into ac- 
count explicitly the scattering anisotropy factor g (with 

J °° P C orr{s)ds = 1): 



Pcorr(s) OC P(s) ■ 1 - f (g) € 



(2) 



mean free path. Both quantities, I and I*, are related by 
I* /I = 1/(1 — <?). The scattering anisotropy parameter g 
is defined as the average of the cosine of the scattering 
angle g = (cos 9). 

To check the validity of Eq. ^ we have performed 
Monte-Carlo simulations of linearly polarized light re- 
flected from a semi-infinite turbid medium (details about 
the simulation method can be found in refs. (0 and l|20|) . 
These simulations use the Mie scattering cross section in 
the range < g < 1 and are able to evaluate numeri- 
cally an exact path length distribution as a function of 
the number of scattering events n and polarization. The 
simulations were done for uncorrelated spherical scatter- 
ers (structure function S(q) = 1) and for a non-reflecting 
interface characterized by z p w I* and z e « 0.7^*. Values 
of Xo/n s = 532 nm for the incident wavelength, n p = 1.59 
for the refractive index of the particle and n s = 1.332 for 
solvent refractive index were used. In Fig. ^ we compare 
the results of the simulations with the prediction of the 
method of images according to Eq. ^ We see clearly in 
this figure that the method of images provides an excel- 
lent description of the path length distribution for the 
case of isotropic scattering (g = 0) , but that the method 
fails to give an equally good description for the case of 



According to Mackintosh and Johni /(l) is of order 
unity and for isotropic scattering (g — > 0) there is no cor- 
rection /(0) = 0. Using the corrected distribution of path 
lengths to fit the simulation results, we have found that 
the function / is well approached by a linear dependancc 
/(<?) = 3g/2 which we assume to be valid also for g < 0. 
As can be seen in figure^ the use of the correction factor 
significantly improves the prediction of Eq. ^ We note 
that, alternatively, for the case of forward-peaked scat- 
tering g ~ 1, other schemes of approximations have been 
suggested. For example the recent Refiii reports that 
the Fokker-Plank equation provides a better description 
than does the (uncorrected) diffusion approximation for 
forward-peaked scattering. 

3. Depolarization length for linear polarization 

Incident polarized light looses its polarization in random 
multiple scatteringtii£i£*i£. For linearly polarized light, 
only two configurations (|| and _L) need to be considered 
(for isotropic samples). Physically, in the || geometry 
more photons are detected for short paths as compared 
to the unpolarized case. In a seminal paper Akkermanns 
et al£ found that the path length distribution for the two 
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configurations can be written as 

P lu (s)=d l± (s)-P(s) 

with the depolarization ratio given by 

7 / \ 1 + 2e- s / e " 
d\\[s) = r„ — 



path length distribution P(s), as given by Eq. has 
been derived for the case of a semi-infinite non-absorbing 

(3) medium. For real systems, both absorption and limited 
container size lead to a loss of photons along a given path. 
For such cases the distribution of path length becomes 
P'(s) — exp(—s/£ a )P(s), where £ a is the characteristic 

(4) absorption length of the medium. In this framework, ab- 
sorption can be taken into account by the modification 



d±(s) 



2 + e" 



(5) 



6t/r -» 6t/r + 3£*/£ a 



(8) 



in terms of the characteristic length of depolarization for 
linearly polarized light £ p . For point like scatterers (g = 
0) Akkermanset al. obtained ^, = £/ln(10/7) ^ 2.804 A 
We find good agreement between EqOHand our numerical 
simulations with £ p as an adjustable parameter (Fig.??). 
For large particles (and therefore large g) the agreement 
is somewhat less good. However polarization effects in 
DWS usually are found weak for g s» 1 and therefore 
we did not attempt to improve the accuracy of Eqs. Q 
and ©(It is worthwhile to note that close to the sam- 
ple surface very interesting polarization effects, such as 
butterfly pattern, persis1*2ii22) . 

In the limit s/£ ^S> 1 Eq. 01 and reduce to : 



1 + 3 - 
- ± — e 

2 4 



■s/e p 



P(s) 



(6) 



We consider this expression the most simple general- 
ization since it captures well intermediate path lengths 
s/i > 3, where polarization effects are important, but at 
the same time the number of scattering events is already 
sufficiently large to apply the diffusion approximation. 

4. Polarization dependence of DWS autocorre- 
lation function 

In transmission geometry the path length distribu- 
tion can be measured experimentally using pulsed laser 
beama^. In reflection however paths are short and there- 
fore the time resolution is usually not sufficient for such 
measurements. An alternative way to probe diffuse light 
propagation is the analysis of temporal fluctuations of the 
scattered light via photon correlation spectroscopy. This 
approach, called diffusing wave spectroscopy (DWS), is 
a sensitive probe to the path length distribution, in par- 
ticular in reflection geometry-'- . The temporal intensity 
correlation function g2{t) which is measured, is related to 
the field autocorrelation function by the Siegert relation 
g\{t) = y/l — g 2 (t) . The latter is directly related to the 
path length distribution: 



gi (t) = / P^e-^f^lt'ds. 



(7) 



In our case the characteristic relaxation time for diffu- 
sive particle motion (diffusion constant D for free Brow- 

nian motion) tq — (k$D) 1 is a known quantity. The 



in Eq. |7J . All our experiments correspond to a case 
where £ a is much larger than £*. 

The solution for the scalar (polarization independent) 
path length distribution is well known 1 (neglecting ab- 
sorption): g x (t) = [e-T^W +e -(7 P +2 7 e M*)] / 2 , where 

x(t) 



\/6£/ T o, 7 P = 
x <C 1 it reduces to 



z p /£* and "f e — z e /£* . In the limit 



91 (t) -- 
0) 



-7 x 



(9) 



7 P + 7e- Note that this 
?* and that for the non- 



with 7 = (z 
expression is independent of 
reflecting interface that we have considered 7 = 1 + 
z e /P* — 5/3. All different paths of length s contribute 
to gi(t). Short paths predominantly contribute to the 
long time decay and long paths contribute to the short 
time decay. Clearly polarization effects modify the path 
length distribution and therefore strongly influence the 
decay of gi(t). Interestingly the shape of gi(t) remains 
more or less unchanged. In most previous studies Eq. El 
though derived for the scalar case, has been applied also 
with polarized light. Although 7 is in principle a well 
defined constant, it has been treated in the literature as 
an adjustable parameter to explain the polarization de- 
pendence of the correlation function. Values of 7j_ 11 in a 
range 1 to 3 have been reported depending on detected 
polarization state, particle size and concentrationiiSiSi. 

Rather than to adjust the parameter 7 we take into 
account the polarization via Eq. I^and find: 



9ix(t) = 



9i\\(t) = 



J P_i(s)e~ 2 ( t / T °) s / r ds 
P±(s)ds 

j p.,( s ) e -2(t/r )s/r ds 

P|| (s) ds ' 



(10) 



Introducing the function h{x) — (e^ lpX ^ + 
)/2, the correlation functions take the 



form (neglecting absorption) 



9i At) 



9i,\\(t) = 



h Mt)] - l_h [ yi (t)] - I h [x 2 (t)] + %h [y 2 (t)] 

h MO)] - fft MO)} - fh [x 2 (0)] + 9 -§h [y 2 (0)} ' 

(11) 

h [ Xl (t)} + \h [ Vl (t)\ - %h [x 2 (t)} - %h [y 2 (t)} 

h MO)] + |A [ yi (0)} - 3 -ih [x 2 (0)] - %h [ W (0)] ' 

(12) 
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d (nm) 


$ (%) 


g 


80 ± 5 


3.9 


0.04 


92 ± 15 


4.1 


0.05 


168 ± 5 


2.0 -> 30.0 


0.29 -> - 


350 ± 12 


1.9 


0.74 


400 ± 8 


1.9 


0.80 


720 ± 14 


1.9 


0.90 


1000 ± 51 


2.0 


0.92 


1500 ± 53 


2.0 


0.92 


114 ± 10 


2.5 -» 7.4 


-0.25 -> - 



0.13 



0.78 



Table 1. Polystyrene spheres used in the experiments. 
Diameter d as obtained from dynamic light scattering, 
volume fraction $ and scattering anisotropy parameter 
g . The latter was obtained either from Mie calculations 
(for hard sphere interactions as described iniS^) or di- 
rect measurements (for the d = 114nm charged spheres 
as described ini&) . 



w here Xi(t) = y/ 6t/T , x 2 (t) = ^6£/t + 3, yi(t) = 
^6t/ To + 3£*/£ p and y 2 (t) = ^6t/r + 3 + 3£*/£ p . 
This set of equations provide a direct relation between 
measurements detecting _L or |j polarized light thereby 
eliminating the need to introduce two adjustable param- 
eters 7_l_ || . We carried out a series of dynamic multiple 
scattering experiments to follow the polarization mem- 
ory of the reflected light intensity. Experiments were 
realized as described m24 but with A = 532nm. The 
sample cells were suspended in a water bath to suppress 
reflections and to maintain a constant temperature of 
T = 22°C. All samples used were made from monodis- 
perse polystyrene particles (n=1.595) suspended in 
water (n= 1.332), except for one case (d — 114nm) where 
we used a mixture of water and ethanoli& (n=1.365). A 
detailed description of all samples is given in tabic ^ 



In Fig. |3 we show a comparison between experi- 
ments and our theoretical expressions. For clarity the 
data is normalized and plotted as a function of x = 
VWto + 31* /la - y/3l*/l a thus removing contributions 
from absorption at x < 0.1 (Eq.©) (The absorption 
length l a has been chosen such that lng\ (x) scales linearly 
at small x- values). We find that the theory describes our 
data very well, with the polarization length £ p being the 
only adjustable parameter. In particular the availability 
of suspensions with negative g values allows a rigorous 
test of the model over the whole interval — 1 < g < 1. 
Note that we use the theoretically predicted values j p , j e 
and do not adjust them to fit the data (as done in all 
previous work). Depolarization lengths obtained from 
fits (Eq. Illfl to the DWS data and from numerical simu- 
lations (Fig.?? ) are presented in Fig. We note again 
that both DWS experiments (<?i_l, <?i ii) are well charac- 
terized by a single £ p even for the most extreme case of 




x=(6t/i o +3l*/l a ) 1/2 - (3I*/I a ) 1 ' 2 

Fig. 3. DWS autocorrelation function for different g - 
values detecting polarized and depolarized light. Lines: 
calculations based on Eq. ^]with £ p adjusted to best fit 
the data, (g = -0.78 : £*/£ a = 0.004, £ p /£* = 5.4), (g = 
0.32 : £*/£ a = 0.0064, l p /l* = 2.21), (g = 0.924, £*/£ a = 
0.016, tp/t* = 0.713). 

In the limit g — we find good agreement with the 
predicted theoretical value of Akkermans £ p /£ = 2.804 5 . 
As the anisotropy parameter g is increased £ p /£ slowly 
increases as well. In the case of forward-peaked scattering 
g ~ 1, we find £ p ~ £* (see also refs. {§) and H)) . This 
means that as g approaches 1 an increasing number of 
scattering events is necessary to depolarize backscattered 
light. Since £ p /£* remains constant for g — > 1 the ratio 
£ p /£ has to increase sharply, as shown in Figure^Jj. In the 
case of backward-peaked scattering g ~ — 1 however, the 
number of scattering events needed to depolarize remains 
virtually unchanged. The characteristic length scale of 
depolarization is still the scattering mean free path £ (and 
not £* !) as in the case of point-like scatterers. Therefore 
it is not surprising to see that in Fig. ^£ p scales as in the 
Rayleigh limit i. e £ p ~ 3£ even for a transport mean free 
path £* ~ £/2. When £ p /£* is plotted as a function of g 
we find an almost linear behavior over the full accessible 
range. At this point our understanding of £ p is limited 
to the particular cases (g = 0, g = 1, and g — —1), a 
more detailed microscopic approach would be required 
to explain the complete dependance of £ p on g, or the 
dependance on other single particle optical properties. 

5. Depolarization ratio of backscattered intensi- 
ties 

An alternative way to study polarization of multiple 
backscattered waves is through the (full intensity) de- 
polarization ratio: 

i\\ + 1± 
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f[P {l (s)+P ± (s)}ds 



P(s) ds , 



from where we identify a s-dependent depolarization ra- 
tio d(s): 



(15) 



Since P||(s) + P ± (s) = P(s) and J °° P(s) ds = 1, there- 
fore: 



d = I P(s)d(s)ds 
lo 



(16) 



5 ( e -lP\fM*/£p + e -(7 P +27e) -\/se*/£ p 

4 V 



In Figure^ we show the experimental values of d. Again 
the agreement is excellent over the full accessible range. 
The data comes very close to the predicted values in the 
three particular cases: d = 0.33 for Rayleigh scattering 
(g = 0), d = 0.14 for forward-peaked scattering (g = 1) 
and d = 0.49 for backward peaked scattering (g = —1). 

6. Conclusion 

In this paper, we have shown how to describe the ef- 
fect of the scattering anisotropy on the depolarization of 
linearly polarized light, and how to use Diffusing Wave 
Spectroscopy (DWS) to determine the characteristic de- 
polarization properties. By means of numerical simula- 
tions, we checked the limit of validity of the diffusion 
approximation when the scattering anisotropy g is in- 
creased, and we have shown how to correct the predic- 
tions by means of an anisotropy dependent cutoff for the 
path length distribution P(s). We discuss for the first 
the time the dependence of the characteristic depolar- 
ization length over the full range of possible values of g 
including the unusual case of negative g values. In our 
description the extrapolation length 7 is a well defined 
constant as required by diffusion theory. Our work thus 
clarifies the meaning of 7, subject to intense discussion in 
the pas1iS^2&. Since our description only uses a single ad- 
justable parameter, it is now possible to fully characterize 
backscattered light with polarization resolved measure- 
ments. We think that this approach can strongly benefit 
applications in the field of soft and bio material analysis, 
as well as diffuse light imaging techniquesA*iii2£. 
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Fig. 4. a) and b) Depolarization length l v from DWS 
measurements and Monte-Carlo simulations. Squares: 
measurements for different particle sizes (random par- 
ticle configuration S(q) = 1). Circles: Strongly inter- 
acting charged particles^. Triangles: Hard sphere data 
from reference 124 Diamonds: Monte-Carlo simulations, 
c) Depolarization ratio directly obtained from the mea- 
sured intensities. Solid lines: calculated from the linear 
fit to l v jl* shown in panel a)^ n 
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